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Abstract 

The Improved Quantum Molecular Dynamics (ImQMD) model incorporated with the statisti¬ 
cal decay model is used to investigate the intermediate energy nucleon-induced reactions. In our 
last work, the description on light complex particle emission has been great improved with a phe¬ 
nomenological mechanism called surface coalescence and emission introduced into ImQMD model. 
In this work, taking account of different specific binding energies and separation energies for various 
light complex particles, the phase space parameters in surface coalescence model are readjusted. 
By using the new phase space parameters set with better physical fundament, the double differ¬ 
ential cross sections of emitted light complex particles are found to be in better agreement with 
experimental data. 
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I. INTRODUCTION 


Since spallation reactions were investigated for the first time 80 years ago by using cos¬ 
mic rays[l], spallation reactions have gained interests in fundamental and applied research 
helds. As optimum neutron source [2], or for energy production and nuclear waste trans¬ 
mutation in accelerator-driven system[3, 4], it recalls the attentions in study on spallation 
reactions. These applications require a large amount of spallation reaction data[5], which 
can not be provided all by experiments [6]. So a theoretical model with powerful prediction 
ability is imperative. The International Atomic Energy Agency (IAEA) and the Abdus 
Salam International Center for Theoretical Physics have recently organized twice interna¬ 
tional workshop to make comparisons of spallation models and codes, including PHITS[7], 
BUU[8], QMD[9], JAM[10], JQMD[11], INCL4[12], ISABEL[13], Bertini[14, 15], Geant[16], 
IQMD[17], RQMD[18], TQMD[19] et al, by merged with various statistical decay models 
such as GEM [20, 21], GEMINI [22], ABLA[23] et ah The productions of the neutron, proton, 
pions and isotopes can be overall described well by the most of given models[5]. But the data 
for the light complex particles (LCPs), i.e. ^H, ^H, ^He, and '^He, can not be reproduced well. 
Only models which have a specihc mechanism to emit energetic clusters, such as coalescence 
during the intranuclear cascade stage or pre-equilibrium emission of composite nuclei, can 
reproduce the high-energy tail of LGPs double differential cross sections (DDXs)[24]. The 
information of yield of hydrogen and helium element is quite important for design of nuclear 
project including target and shield. Therefore, the description of model on LGPs emission 
should be improved to satisfy the request of the applications[24, 25]. 

In our previous works, by applying the Improved Quantum Molecular Dynamics 
(ImQMD) Model merged statistical decay model, a series of studies on the proton-induced 
spallation reactions at intermediate energies have been made [26-28]. Nuclear data includ¬ 
ing neutron DDXs, proton DDXs, mass, charge and isotope distributions can be overall 
reproduced quite well. In our last work[29], with a phenomenological mechanism called sur¬ 
face coalescence and emission is introduced into ImQMD model, the description on LGPs 
produced in nucleon-induced reactions are great improved but still not good enough. The 
motivation of this work is to further improve the surface coalescence and emission model. We 
expect the model with new parameters can describe the experimental data more accurately. 

The paper is organized as follows. In the next section, we make a brief introduction on 
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the model. In section III, we give the calculation results and make some discussion. Finally, 
a summary is given in Sec. IV. 

II. MECHANISM OF SURFACE COALESCENCE 

ImQMD model merged GEM2[20, 21], the same method as one used in our last work[29], 
is used in this work. The very detailed introduction on the ImQMD model can be found 
in ref. [29] and the reference therein. Here we only make an introduction on the surface 
coalescence mechanism which has been introduced into ImQMDOd model. After the incident 
nucleon has touched the target nuclei to form the compound nuclei, we define a sphere core 
with radius Rq surrounding by a surface with width Dq. At each time step, any fast nucleon 
passing the surface region to leave the compound system is taken as leading nucleon. An 
inspection is made over all other regular nucleons in order to check whether there are one or 
several nucleons close enough in phase space to allow the formation of a stable composite- 
particle satisfying the following condition 

Rim ^ ^im ^ ^0 witll Rim ^ f fm, 

Where the Rim and Pim are the Jacobian coordinates of the Ah nucleon, i.e., the relative spa¬ 
tial and momentum coordinates of the considered Ah nucleon with respect to the subgroup 
cluster m. The value of ho is adjustable by htting the experimental data. And the condition 
Rim > 1 fai gets rid of the unreal LCPs constituted by the nucleons being too close to each 
other, due to the repulsion of nucleon-nucleon interaction in short distance [30]. In present 
work, the following LCPs are considered: ^H, ^H, ^He, and ^He. By the method of LCPs to 
be constructed, candidate nucleon belongs to a heavier LCP also belongs to a lighter LCP. 
So LCPs are checked to be emitted according to the priority list: '^He>^He>^H>^H, say the 
heavier LCP is first tested for emission, as the same order as that in Ref. [30-32]. Finally the 
candidate LCP can be emitted or not depends on whether its kinetic energy is high enough 
to tunnel through the Coulomb barrier. The kinetic energy of the LCPs can be calculated 
as 

-^Icp 

-E'lcp + Vi) + 

Bicp ( 2 ) 

i=l 

Where the Ei and V are the kinetic energy and the potential energy of the Ah constituent 
nucleon, respectively, Ai^p and are the mass number and the binding energy of the 
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LCP, respectively. If all conditions are meet, the LCP is emitted in the direction of its c.m. 
momentum. Otherwise, all nucleons in the “LCP” are set free and become available again in 
the nucleus and in the ImQMD process, the leading nucleon is emitted as a free nucleon. In 
this coalescence model, composite-particles are thus not allowed being formed in the interior 
of the nucleus but only in the surface layer. It is reasonable according to the knowledge from 
nuclear structure and reaction. For each leading nucleon, the LCP formation and emission 
are tested in the priority list. At each time step, the same test is repeated until the end of 
ImQMD simulation. 

III. RESULTS AND DISCUSSION 

By systematic comparison between calculation results and experimental data of nucleon- 
induced reactions, the parameters in the surface coalescence model are hxed. One best 
choice is Ro= 1.4A^/^, Dq=2.3 fm, and phase space parameters 

{ 200 fm MeV/c, A'lab < 300 MeV, 

260 fm MeV/c, 300 MeV < i?iab < 500 MeV, (3) 

330 fm MeV/c, Ei^h > 500 MeV. 

We call this phase space parameters set as setl. Then with the hxed parameters, chosen 
once for all, the prediction power of the model is tested by the nucleon-induced reactions 
on various targets with energies from 62 to 1200 MeV. And it is found that, with surface 
coalescence mechanism introduced into ImQMD model, the description on the DDXs of LCPs 
is great improved. The experimental data can be overall reproduced well. But the calculation 
results are still a little far from satisfactory in detail. The experimental data of each single 
LCP emitted in the same reaction can not be simultaneously and satisfactorily reproduced 
with one given hg. For example, calculated DDXs of and in the reactions n-l-^^Cu at 
317, 383, 477 and 542 MeV are represented in Fig. 1, with ho=200, 260 and 330 fm MeV/c 
adopted, respectively. Calculation results without surface coalescence mechanism are also 
shown in the hgure. As a rule in this article and for the sake of clarity, DDXs are displayed 
after multiplication by 10°, 10“^, 10“^, 10“°, etc., as noted in the hgures. One can see that, 
for ^H, larger ho is required to reproduce the experimental data. But for °H, calculations 
with large ho overestimates the experimental data. Finally, the compromised values of ho 
listed in Eq. (3) are selected to make calculations to overall reproduce experimental data. So 
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FIG. 1: (Color online) Calculated DDXs of and in the reactions n+®^Cu at 317, 383, 477 
and 542 MeV. /io=200, 260 and 330 fm MeV/c are adopted, respectively. Experimental data are 
taken from [33]. 

it seems that, to reproduce the experimental data of each LCPs as well as possible, various 
ho should be adopted when construct different LCPs. 

From the knowledge of nuclear structure, we know that the nuclei with higher specific 
binding energy is more stable. In table I, the specific binding energies for LCPs are repre¬ 
sented. The descending order of specific binding energies for LCPs is ^He>^H>^He>^H. It 
means that the nucleons in ^He are bound most tightly and the nucleons in are bound 
most loosely. On the other hand, separation energy also reflects the tightness between the 
nucleons in LCP. From Table I, one can see that the descending order of separation energies 
for LCPs is ^He>^H>^He>^H. It means that, compared to nucleon in ^H, the last nucleon in 
^He is more difficult to be separated from mother nuclei. In surface coalescence model, one 
by one candidate nucleon is coalesced to leading nucleon to form LCP, so the condition Eq. 
(1) for the last candidate nucleon should be adjusted according to the effect of specific bind¬ 
ing energies and separation energies of LCPs. The LCP with large specific binding energy 
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TABLE 1: Specific binding energies, separation energies and value of /iq for LCPs, respectively. 


LCP 

construction 

Eb/A (MeV) 

Sn,p (MeV) 

/io(fm MeV/c) 

2H 

p + n 

1.1 

2.24 

300 


n-|-^H 

2.9 

6.24 

260 

^He 


2.6 

5.48 

230 

^He 


7.1 

19.82 

200 

^He 

n-|-^He 

7.1 

20.58 

200 


and separation energy requires the last candidate nucleon to be more close to the subgroup 
cluster in the phase space with smaller h^. With Rq= Dq=2.3 fm, as same as ones 

used in last work, the value of ho are determined by htting the DDXs experimental data of 
^H, ^H, ^He, “^He in the nucleon-induced reactions on various targets with incident energies 
from 62 MeV to 1.2 GeV. In table I, we list the value of ho used to construct different LCPs. 
For example, a leading proton can pick any neutron satisfying Rim x Pim < 300 fm MeV/c 
to form ^H. But for a generated ^He to form ^He, it must pick up the last neutron satisfying 
Rim X Pim < 200 fm MeV/c. According to > S'p(^He), ho(^H) should be smaller than 

ho(^He). But one can see that ho(^H) > ho(^He) is taken, because there is only difference of 
0.76 MeV in separation energies between two particles, and the emission priority >3He>3H 
is also a factor to be considered. The similar phase space parameter set are used in Ref. 
[30]. We call this phase space parameters set as setll. 

Figures 2-4 represent calculated DDXs of LCPs in several nucleon-reduced reactions with 
parameter setl and setll adopted, respectively. One can see that, although for some LCPs, 
for example, in 175 MeV p-|-^®Ni and ^He in 1200 MeV p-|-^®^Au reaction, the results 
calculated with parameter setll are not as good as the results calculated with parameter 
setl, the most experimental data can be described better with parameter setll. Because by 
using parameter setll, the independent phase space parameters for each LCP is helpful to 
rather satisfactorily £t the experimental data. 

Finally, reactions of 317, 383, 425, 477, 542 MeV neutrons hit ^°®Bi are used to test 
the model. The comparisons between calculated DDXs of LCPs and experimental data are 
illustrated in Fig. 5. From the hgure, one can see that excellent agreement of LCPs DDXs 
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FIG. 2: (Color online) Calculated DDXs of and in the reactions n+®^Cu at 317, 383, 425, 
477 and 542 MeV with parameter setl and setll adopted, respectively. Experimental data are taken 
from [33]. 




FIG. 3: (Color online) Calculated DDXs of LCPs in the reactions n+^®Fe at 175 MeV (experimental 
data are taken from [34]), and p+^®Ni at 175, 1200 MeV (experimental data are taken from [35, 36]) 
with parameter setl and setll adopted, respectively. 
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FIG. 4: (Color online) Calculated DDXs for LCPs in the reactions p+^^Al, p+^^Fe (middle pane) 
at 62 MeV, (experimental data are taken from [37]), andp+^®^Au at 1200 MeV (experimental data 
are taken from [38]) with parameter setl and setll adopted, respectively. 


are achieved. 


IV. SUMMARY 

In our last work, a phenomenological surface coalescence mechanism is introduced into 
ImQMDOS model. With surface coalescence mechanism introduced, the description on the 
light complex particles emission in nucleon-induced reactions is great improved, the exper¬ 
imental data of double differential cross sections of LCPs can be reproduced well. In this 
work, the phase space parameters in the surface coalescence model are readjusted with spe- 
cihc binding energy and separation energy effects for various LCPs being considered. The 
experimental data are better accounted with new parameters set. And what more important 
is that the new parameters set has better physical fundament. 
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FIG. 5: (Color online) Calculated DDXs of light complex particles produced in the reaction n+^^®Bi 
and re+^^^Bi at 317, 383, 425, 477 and 542 MeV, respectively. Experimental data are taken from 
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